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Abstract 

In order to estimate the possibility to observe exotic physics in a neutrino telescope, it is 
essential to first understand the flux of atmospheric neutrinos, muons and dimuons. We 
study the production of these leptons by high-energy cosmic rays. We identify three main 
sources of muons of energy E > 10 6 GeV: the weak decay of charm and bottom mesons and 
the electromagnetic decay of unflavored mesons. Contrary to the standard assumption, we 
find that 77 mesons, not the prompt decay of charm hadrons, are the dominant source of 
atmospheric muons at these energies. We show that, as a consequence, the ratio between 
the neutrino and muon fluxes is significantly reduced. For dimuons, which may be a 
background for long-lived staus produced near a neutrino telescope, we find that pairs of 
E ?z 10 7 GeV forming an angle above 10~ 6 rad are produced through D (80%) or B (10%) 
meson decay and through Drell-Yan proceses (10%). The frequency of all these processes 
has been evaluated using the jet code PYTHIA. 
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1 Introduction 



We observe a flux of cosmic rays (protons free or bound in nuclei) that extends up to 
energies of 10 11 GeV [1]. When a very energetic cosmic ray enters the atmosphere it will 
collide with a nucleus of air at an approximate altitude of 20 km. The interaction breaks 
the primary proton, starting an air shower that includes millions of particles, collisions 
and decays. The development of the shower along the atmosphere (number of particles 
and energy deposited at different depths) is fairly well understood [2,3]. In particular, 
computer simulations provided by codes like AIRES [4] or CORSIKA [5] reproduce well 
the profile of the shower and the number of charged particles (mostly muons) that reach 
the ground. 

Most of the processes that take place inside the air shower are soft hadronic collisions 
and particle decays. Typically, a hadron is initially produced at a given point in the 
atmosphere with a large boost towards the ground. If this hadron is a proton or a neutron 
it will move along a hadronic interaction length (Ai n t ~ 1-6 km, depending on the altitude) 
and will collide with a nucleus, exchanging momentum with q 2 < 1 GeV 2 . The process 
will result into a leading hadron plus several other hadrons sharing the total energy. On 
the other hand, if the initial hadron is a charged pion (or a kaon) the sequence of events 
depends critically on its energy E. For E larger than E^ lit = 115 GeV [6] (or E^. it = 855 
GeV) its interaction length is smaller than its decay length, the pion tends to collide 
and behaves much like the stable proton. Q For smaller energies, however, it decays and 
produces the conventional muon and neutrino fluxes. Finally, a third possibility is that the 
initial hadron decays into other hadrons, photons or charged particles very fast through 
strong or electromagnetic interactions (this is the case for resonances, neutral pions, and 
other hadrons that do not decay weakly). 

Some physical observables, however, may depend on rare processes of higher q 2 that also 
occur inside the shower and cannot be overlooked. For example, consider the atmospheric 
muon flux above 10 6 GeV. At these energies charged pions collide with the air before they 
can decay. The origin of these muons, as it has been widely described in the literature 
[7-14], is the prompt decay of charmed hadrons produced in hadron- nucleon collisions. 
Charm production involves q 2 > (2m c ) 2 and is less frequent than pion or kaon production. 
One may argue that the number of these very energetic muons is negligible compared to the 
total number of muons reaching the ground, and also that the energy that they take from 
the shower is a non-significant fraction of the total energy deposited in the atmosphere. 
Such arguments could justify the absence of charm and bottom hadrons in the air-shower 
simulations performed by AIRES or CORSIKA. Even if negligible there, however, these 
muons may be observable in a neutrino telescope like IceCube [15,16], specially from near- 
horizontal directions. Notice that given a zenith angle 9 fixing the ice column density 
faced by an atmospheric particle in its way to the telescope, only muons above some 
energy threshold can reach the detector. An IceCube measurement of the muon flux at 
E > 10 6 GeV would provide information about hadronic collisions which is complementary 

TThe actual value of E^ Tit depends on the air density at the point where the pion is produced. 
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to collider data, as this forward physics at such high energies is not available there. 

In this paper we evaluate the production of high-energy atmospheric leptons: muons, 
neutrinos and dimuons. We are interested in the energy region around 10 GeV, where a 
D meson tends to interact before decaying just like a pion does already at lower energies. 
Our interest is motivated by several factors. The atmospheric muon and neutrino fluxes 
are correlated, so a measurement of the former determines the latter. Ultrahigh energy 
atmospheric neutrinos are of interest because (i) they are a background to a possible cosmic 
neutrino flux to be observed in IceCube, and (ii) it has been suggested [17] that they could 
be a source of new physics {e.g., uN — > ffX) observable at a neutrino telescope. Moreover, 
atmospheric muon pairs from quasi-horizontal directions are themselves a background to 
exotic new physics at a neutrino telescope. For example, they could be confused with a pair 
of long-lived charged massive particles (CHAMPs) produced near the detector through a 
vN interaction [18]. 

We think that previous estimates of the muon and neutrino fluxes at ultrahigh energies 
have overlooked some very relevant effects. Most notably, they do not include muons 
produced in the electromagnetic decay of unflavored mesons. In addition, some of them 
do not properly include the propagation in the atmosphere of charmed hadrons of E > 10 7 
GeV (see below). Finally, we include in our study the Drell-Yan process qq — > 7/Z — > /j,fj,, 
which cannot be neglected when studying the background to exotic physics at neutrino 
telescopes. 



2 Atmospheric muon and neutrino fluxes 

We will focus on leptons of energy above 10 6 GeV. As explained before, at these energies the 
contribution from charged-pion decay is negligible. Actually, we will take secondary pions 
and kaons at these energies as stable particles that may produce muons in their collision 
(not decay) with a nucleus of air. In [19] we find an estimate of the total atmospheric 
flux of primary plus secondary nucleons (N) and long-lived mesons (n ± and K). We have 
used the jet code PYTHIA [20] to simulate their collisions and identify the main sources 
of leptons, which are the following. 

(i) The standard source of muons of energy above 10 6 GeV is the prompt decay of charmed 
hadrons. These hadrons decay into muons with a branching ratio of about 10%. 

(ii) As the energy of a charmed meson grows, its decay length Ad ec becomes larger than 
its interaction length Ai nt , reducing the probability of decay into leptons before colliding. 
The possibility that the charmed hadron interacts several times in the air and then decays, 
however, can not be neglected. The main reason is that a D meson does not propagate in 
the air like a pion. The charmed hadron is much more penetrating. We estimate that in a 
hadronic interaction it loses just a fraction 

x = A~o.3 (2.1) 

m c 
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of the energy lost by a pion of equal energy. After each interaction with air there is 
always a leading charmed hadron carrying more than 70% of the initial energy. We find, 
for example, that the decay of charmed hadrons that have interacted with the atmosphere 
before increases the prompt lepton flux at E > 10 GeV by 30%. Previous analyses [9] have 
assimilated the propagation of a D meson (or a A c baryon) to that of a pion (or a A), which 
reduces the energy of the leading hadron and makes this contribution negligible. Other 
studies however, have included this effect properly (for example, through a regeneration 
factor Z cZ = 0.8 in [21]). 

(in) The two processes above would also work for bottom hadrons, which have similar 
lifetime and semileptonic branching ratio. Of course, B mesons are less frequently produced 
than D mesons by cosmic rays in the atmosphere, but their decay length is smaller. Their 
contribution to the 10 7 GeV muon flux is a 10% of the one from charm hadron decays. 

(iv) The electromagnetic decay of unflavored hadrons, most notably rj mesons, has been 
almost completely overlooked in the literature (the only generic mention that we have found 
is given by Ryazhskaya, Volkova and Zatsepin in [22]). Since these hadrons do not contain 
any heavy quarks, they are more abundant in air showers than charmed or bottom hadrons 
(rj mesons are produced through low q 2 interactions, just like pions or kaons). In addition, 
since they decay through electromagnetic processes their lifetime is much shorter than the 
one of charged pions, kaons or D mesons, whose decay is mediated by weak interactions. 
The decay length of a 10 10 GeV rj meson, for example, is just 280 m. On the other hand, 
the preferred decay modes are into photons and neutral pions, the branching fraction to 
yU + /i _ 7 (around 3.1 x 10~ 4 ) being suppressed by an cxqed factor. In the lenguage of the 
Z-moment method, these production and decay rates would translate into an approximate 
Z pr] Z p k l and Z VI1 ~ 2 x 10~ 3 Zk l11 - Notice also that, in contrast with the other (weak) 
processes, these contributions break the correlation between the atmospheric muon and 
neutrino fluxes. 

(v) Finally, we will also include in our analysis the Drell-Yan process qq — > j/Z — > 

This scattering will be relevant to determine the flux of muon pairs (dimuons) forming a 
minimum angle such that they can be resolved at a neutrino telescope. 

Our results are summarized in Figs. 1-2. We have simulated hadron (h = p,n,TT + ,K 
and antiparticles) collisions with atmospheric nucleons and identified the muons and neu- 
trinos produced through the five processes described above. We have estimated the prob- 
ability for a process X by comparing its cross section (obtained from PYTHIA) with 
the total cross section with the air: 

A n hN 

-P h x(E) - , (2.2) 

where A = 14.6 is the averaged atomic mass of a nucleus of airH PYTHIA has simulated the 

* By using the factor A instead of A 2 / 3 we ignore the screening between target nucleons and include 
the possibility that the incident hadron collides with more than one nucleon inside the nucleus of air (its 
partons do not disappear after the first interaction). As a consequence, (|2.2[) expresses the number of 
(short distance) interactions per collision with a nucleus, and could be larger than one. 
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Figure 1: Muon flux started by nucleon and meson collisions from the sources (i-iv) 
(solid) and just from meson collisions (dashes). Conventional muons from pion decay are 
not included. 



production of five types of charmed hadrons (D°, D + , D+, A+, and their antiparticles), 
seven of bottom hadrons (B°, B + , Bf, B£, A°, S^~) and three of unflavored mesons 
(rj, p and <fi). PYTHIA also takes care of gluon emission, the decay of hadronic resonances, 
and the decay into leptons of the parent heavy or unflavored hadron. These probabilities 
have then been convoluted with the total fluxes in [19]. 

We have estimated for all these hadrons an interaction length in the air of 5 km, 
allowing that they decay between interactions (we have neglected the variation of Xd ec with 
the altitud; since most of these long-lived hadrons are produced between the second the 
fourth interaction lengths, the 5 km should be an acceptable approximation). We have 
assumed that in each interaction they lose 30% (charm) or 10% (bottom) of their energy. 

In Fig. 1 we plot the flux of muons at different energies. We have separated the 
contributions from the prompt decay of charmed hadrons (that may have interacted before 
decaying) (i-ii), from B decays (Hi), and from 77 decays (iv). We have also separated 
the contribution from collisions of secondary mesons (charged pions and kaons) with air 
nucleons, which is about a 15% of the total. 

One important point concerns the procedure used to estimate the probability paec that, 
once produced, a hadron decays before interacting. PYTHIA assigns to the hadron a 
proper lifetime r that is distributed around the average value To. We boost r and calculate 
the probability that in this time the charmed hadron does not interact with the atmosphere 
before decaying: 

Pdcc = exp <^ — \ , (2.3) 

with Ai n t ~ 5 km. Notice that always taking the average value r instead of r, the decay 
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E e [GeV] E v [GeV] 

Figure 2: Total muon and neutrino fluxes (left) and relative contribution from charmed 
and bottom hadrons in the neutrino fluxes (right). 

probability of a charmed meson of energy above E cr - lt = mX mt /ro would be much smaller, 
since one would exclude the possibility of r being significantly smaller than tq. We think 
this may be the reason why other estimates predict a sharper drop in the muon flux at 
energies above 10 8 GeV. Other than that, the contribution from charm production that 
we obtain with PYTHIA (based on perturbative QCD) is similar to the one given in [9]. 
Soft physics is treated by PYTHIA using reggeon and pomeron exchange [23, 24] and the 
string fragmentation model [25]. 

In Fig. 2 (left) we show the total muon flux along with the fluxes of the three different 
neutrino species. The ratio of muon neutrinos to muons changes from 1 to 0.2 at E = 10 8 
GeV with the inclusion of the electromagnetic decays of unflavored hadrons. In Fig. 2 
(right) we give the contribution of charmed and bottom decays to the three neutrino 
species. The contribution to the tau neutrino flux from the prompt decay of D s mesons 
that we obtain is smaller than what one would expect from a simple estimate (notice that 
other charmed hadrons are too light to decay into tau leptons) . This result does not agree 
with previous analyses [21,26-28], and seems to be caused by a too low decay rate of D s 
into tv t given by our PYTHIA simulation! 

3 Atmospheric flux of muon pairs 

Two muons separated by a distance of 50 meters crossing a neutrino telescope could look 
similar to a pair of CHAMPs. If the CHAMPs come upwards [29] , their origin could only be 
a neutrino interaction at some distance from the detector: the upward atmospheric dimuon 

§ We have used the default PYTHIA 6.418 (fortan) distribution of the code. 
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Figure 3: Total dimuon flux (left) and dimuon flux requiring a minimum angle of 10~ 6 rad 
(right). We separate the contribution from charm decay, bottom decay, unflavored meson 
decay, and Drell-Yan processes. 

background vanishes whereas the typical distance (angle) between two muons produced 
by a neutrino near the detector is always smaller. However, if the CHAMPs reach the 
telescope from a quasi-horizontal direction {6 = 70°-95°), they could have been produced 
in the atmosphere (in a hadron-nucleon interaction) or in the Earth (a neutrino-nucleon 
collision), and the evaluation of the expected number of atmospheric muon pairs becomes 
essential. 

We have estimated this flux at very high energies E^. We require that the total energy 
is balanced, with each muon carrying at least 1% of E^. This is necessary as both muons 
must be able to reach the telescope from a given distance. Their origin are also the five 
processes discussed in the previous section. Heavy quarks are always produced in pairs 
by cosmic rays, so it may be that both D (or B) mesons decay into muons and define a 
dimuon. The decay of an 77 meson gives dimuons (never single muons), but the transverse 
momentum is very small and does not provide enough separation to resolve them. In 
addition, the Drell-Yan process becomes important when a minimum angle between the 
two muons is imposed (notice that 77 production and decay into muons, qq — > 77 — > 
can be interpreted clS db Drell-Yan at q 2 ~ ml). 

In Fig. 3a-b we plot the total dimuon flux (left) and the flux requiring a minimum 
angle of 10~ 6 rad (right). We separate the contribution from charm decay, bottom decay, 
unflavored meson decay, and Drell-Yan processes (which include 7 and Z exchange). It is 
apparent that the requirement of a minimum angle cuts all dimuons from 77 decays, while 
prompt leptons dominate up to 10 7 GeV and Drell-Yan at higher energies. In Fig. 4 we 
give the total flux after cuts of 0, 10~ 8 , 10~ 6 and 10~ 4 rad in the angle between the two 
muons. 
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Figure 4: Total dimuon flux for different cuts in their opening angle. 

4 Summary and discussion 

Cosmic rays produce extensive air showers with millions of particles of different energy. 
The atmosphere acts as a calorimeter and absorbs most of the initial energy. However, a 
small fraction of this energy goes to very penetrating particles that are able to enter the 
ground and reach a neutrino telescope like IceCube: muons and neutrinos. If the detectors 
there look down, a signal could only be due to neutrino interactions near (or inside) the 
telescope. The interacting neutrino could be atmospheric or cosmogenic, whereas the signal 
could correspond to standard (muons crossing the telescope upwards) or exotic (stau pairs 
from the same direction). In any case, the determination of the origin of the signal requires 
an accurate estimate of the atmospheric neutrino flux at these energies [30] . 

In addition, if the telescope also covers quasi-horizontal zenith angles (8 = 70°-95°) 
there appear several new and interesting possibilities that could explain the origin of a 
signal. The large ice column density along these directions filters most of the components 
in an air shower, just like from upward directions. Actually, everything but neutrinos and 
very energetic muons will be absorbed by the ground. The experimental determination 
of the muon flux from these inclinations would be a direct test for the different computer 
codes used to simulate air showers, and also an indirect measurement of the atmospheric 
neutrino flux at these energies. More importantly, these horizontal directions are the most 
promising in the search for exotic physics, since this physics may have been produced 
not only by cosmogenic or atmospheric neutrinos in the ice, but also in muon-nucleon or 
hadron-nucleon collisions in the atmosphere. For example, a couple of massive, long-lived 
(very penetrating) staus produced in the air via gluino decay would easily reach the center 
of a neutrino telescope from horizontal directions, but not from below, as they would be 
unable to cross the Earth. The identification of this type of exotic physics requires an 
accurate estimate of the muon and dimuon atmospheric background. 
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We have determined the flux of muons and neutrinos at energies above 10 6 GeV. We 
have found that r\ decay, not the prompt decay of D mesons, is the main source of muons at 
these energies. Our results confirm the estimate given in [22]. D and B mesons decay via 
W exchange and become long-lived at very high energies just like charged pions do at lower 
energies. In contrast, 77 mesons decay fast through electromagnetic interactions (just like 
neutral pions). The rj branching ratio to muons is suppressed by an oqed factor, but this is 
compensated by the fact that they are more abundantly produced in air showers than charm 
or bottom hadrons. An important consequence of these unflavored meson contributions 
is that the atmospheric neutrino to muon ratio is smaller than the one obtained by other 
authors. 

We have also determined the flux of very energetic dimuons. These muon pairs would 
be a background to long-lived CHAMPs produced in the air or the ground and reaching a 
neutrino telescope from horizontal directions. If a minimum angle between the two muons 
(so that they can be resolved) is imposed, we have obtained that the dominant sources are 
heavy hadron decays and Drell-Yan processes. 

We think that neutrino telescopes may provide opportunities that go beyond the study 
of neutrino interactions at ultrahigh energies. The study of the muon flux from near- 
horizontal directions would provide important information about hadronic collisions at 
energies and q 2 (forward direction) not accessible in colliders, and it could even reveal 
signals of non-standard physics. 
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